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Abstract: The discovery of a Higgs boson with a mass of 126 GeV at the LHC when 
combined with the non-observation of new physics both in direct and indirect searches 
imposes strong constraints on supersymmetric models and in particular on the top 
squark sector. The experiments for direct detection of dark matter have provided with 
yet more constraints on the neutralino LSP mass and its interactions. After imposing 
limits from the Higgs, flavour and dark matter sectors, we examine the feasibility for 
a light stop in the context of the pMSSM, in light of current results for stop and other 
SUSY searches at the LHC. We only require that the neutralino dark matter explains 
a fraction of the cosmologically measured dark matter abundance. We hnd that a 
stop with mass below ~ 500 GeV is still allowed. We further study various probes of 
the light stop scenario that could be performed at the LHC Run - H either through 
direct searches for the light and heavy stop, or SUSY searches not currently available 
in simplihed model results. Moreover we study the characteristics of heavy Higgs for 
the points in the parameter space allowed by all the available constraints and illustrate 
the region with large cross sections to fermionic or electroweakino channels. Finally we 
show that nearly all scenarios with a small stop—LSP mass difference will be tested by 
XenonlT provided the NLSP is a chargino, thus probing a region hard to access at the 
LHC. 
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1 Introduction 

The Large Hadron Collider (LHC) has completed its first run (Run-I) with an unprece¬ 
dented success. A Higgs particle has been discovered with a mass ~ 126 GeV [1-6]. 
Its couplings to the Standard Model (SM) electroweak gauge bosons have been estab¬ 
lished to be close to the SM expectations by measurements of signal rates [6-15] and 
spin-parity determinations [7, 16-18]. While there remains considerable room for devi¬ 
ations in the couplings to fermions, no sign of any New Physics (NP) has been detected 
yet. The search for new states near the electroweak scale has also been frustrating. 
The trudge of null results has shrunk the parameter spaces of weak scale NP models 
considerably, and models of weak scale Supersymmetry (SUSY) are not exceptions. 

Indeed, the general tone at the moment is fairly lugubrious for SUSY enthusiasts. 
The simplest versions of the Minimal Supersymmetric Standard Model (MSSM), with 
the simplest assumptions about the high-scale theory, are under increasing tension with 
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a wide range of experimental data inclnding the ~ 126 GeV mass of the Higgs boson [19- 
22] . However, the above argnment can be tnrned aronnd to advocate many convincing 
reasons to stndy low energy SUSY. For example, that the electroweak symmetry is 
broken by an elementary scalar whole mass is below 135 GeV, is, in fact, a prediction 
of the MSSM. One shonld also be reminded that the MSSM has excellent deconpling 
properties which keep electroweak precision observables nnder control and to some 
extent also ameliorate the tension with limits on flavonr observables. Moreover, the 
states which are directly related to the natnralness of Electro-Weak Symmetry Breaking 
(EWSB) are not constrained severely by the direct searches yet. 

Existence of these states aronnd the TeV scale is implied by demands of natnralness, 
which is a mnch discnssed issne in the context of Snpersymmetric theories (for example 
see [23, 24]). The main point can be nnderstood by considering the issne of stabilization 
of the Higgs mass against radiative corrections. To be specihc, the correction to the 
Higgs mass dne to radiative effects can be written as, 

*nJ(AEw)~Ai„svl.i(^) , (1.1) 

V ^'■EW / 

where A^ess denotes the scale at which SUSY breaking effects are mediated to the 
MSSM and a common mass scale Asusy for all the SUSY particles has been assnmed. 
Eq. (1.1) immediately makes it clear why SUSY particles (especially those which conple 
strongly to the Higgs) < TeV are desired. 

For moderate to large tan(3 = {Hu)/{Ha), e.g. tan/? > 2, the Higgs mass in the 
MSSM can be written as [23] 


'ml — —2(|/ip -I- mjjjtree + 'mjjjrud), 


( 1 . 2 ) 


where fi is the snpersymmetric Higgs mass parameter, and and are 

the tree-level and radiative contribntions to the soft SUSY breaking mass sqnared for 
Hu- The dominant radiative correction to proportional to the top qnark Ynkawa 
conpling is given by [23], 


m 


H,, 


M 
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(1.3) 


where yt is the top Ynkawa coupling, m~ and m~ are the soft SUSY breaking mass 
sqnared parameters for the third-generation sqnark donblet and singlet np-type sqnark. 
At is the scalar trilinear interaction parameter for the top sqnarks ^ and My denotes 

^Note that we use stop and top squark interchangeably. 
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an average mass scale for the top squarks. Recall that in Supersummetric theories a 
light Higgs is ‘natural’ in the sense that the stabilization of the Higgs mass around 
the EW scale is guaranteed by the symmetry. The destabilizing effects come from 
SUSY breaking. If none of the terms on the right-hand-side of Eq. 1.2 are much larger 
than the left-hand-side then it implies that no hne tuning of parameters in the theory 
is needed to guarantee the low Higgs mass. Thus the amount of ‘cancellations’ (hne 
tuning) required to satisfy Eq. 1.2 is then a measure of ‘naturalness’. As an example if 
we dehne the hne-tuning parameter A = 2m|^^ |rad/'^h upper limit on the 

top squark mass scale as a function of the hne-tuning parameter [25-27], 



(1.4) 


Upper bounds on /i and the gluino mass can also be obtained in a similar way, 



(1.5) 


( 1 . 6 ) 


Equation 1.6 follows from taking into account, in the leading-logarithm (LL) approx¬ 
imation, corrections to m\ coming from gluino mass. These come from the gluino 
induced corrections to the Higgs potential arising at two loops which in turn come 
from gluino corrections to the stop mass at one loop. 

If we forget about the Higgs mass and other indirect constraints for a moment, then 
from direct searches alone one can still have ~ ~ 500 GeV and ~ 0. This 

amounts to a tuning A < 10 (assuming Amess = 20 TeV). However, the direct search 
bound on the gluino mass m-g > 1500 GeV requires, according to Eq. 1.6, the tuning 
to be A > 15. This means that, as far as the direct search bounds are concerned, it is 
the gluino mass that has stronger effect on hne tuning than the stop mass. Of course, 
when constraints from the Higgs mass are taken into account, low values for At are not 
allowed (if at least one of the top squarks is desired to be light). This makes the tuning 
much worse (A > 50). 

The effect of gluino mass on hne-tuning can also be understood if one considers the 
running of the stop mass. Indeed the top squark is a scalar and its mass is subject to the 
same hne tuning problem as the Higgs mass. More precisely, the leading contribution 
is given by [23] 


dini^/ij iOTT” o 



(1.7) 
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This means that the stop mass is attracted towards the gluino mass at low energies and 
thus a light stop does not seem very ‘natural’ in view of the rather high lower bounds 
on the gluino mass implied by LHC Run-I data. One way to allow a light stop and 
still be consistent with naturalness, is to somehow weaken the rather strong bounds on 
the gluino mass implied by the current LHC data. This can be done in SUSY models 
like compressed SUSY, stealth SUSY and R-parity violating SUSY, hence the renewed 
interest in such models [28-33]. 

One must further appreciate that the naturalness criteria, to some extent, tend to 
be subjective. In fact. Ref. [34], argues that an appropriate model-independent measure 
of hne tuning is not in terms of the high scale quantities, but rather in terms of /i. Thus 
a light Higgsino is the one robust demand one can make on the particle spectrum by 
requiring naturalness, without making any assumptions about the high scale physics. 
In this analysis they construct a measure of naturalness which is ‘independent’ of 
the precise model of the high scale physics. This measure can be small even when 
the conventional measures of hne tuning such as A take large values. Due to such 
subjectivity in the naturalness criteria, in our analysis we will cover values of parameters 
beyond the nominal upper bounds indicated by Eqs. 1.4-1.6, for a given value of A. 

Finally, irrespective of the amount of hne tuning introduced it is interesting to 
investigate to which extent light stop masses are allowed by all the existing direct and 
indirect measurements, and to see how well this region can be tested at the LHC. Here, 
we will work within the framework of the CP and havour conserving phenomenological 
MSSM (pMSSM) with 19 free parameters [35, 36]. Within this model we will concen¬ 
trate on the electroweak scale parameters that are most relevant for the Higgs and dark 
matter (DM) sector, that is the gaugino, Higgs and third generation squark parame¬ 
ters. We then explore the parameter space of the model allowed by havour constraints, 
Higgs properties and Higgs searches, direct LHC searches for SUSY particles as well as 
DM constraints (relic density and direct detection). 

The collider, havour and DM constraints on the general pMSSM were explored in 
several publications [37-45] and the impact of stop searches was also considered [41, 
43, 44, 46-56]. Possible probes to improve the bounds on stops and sbottoms with 
new observables and/or better background reduction have been explored in several 
works [57-73] and the difficult case of a compressed spectrum has been considered [74- 
84] . In this work we concentrate specihcally on the parameter space of the model where 
a light stop is allowed and we rely on simplihed model constraints on the stop mass 
at LHC Run-I. The aim is to understand the constraining power of the current search 
results and demonstrate possible ways in which these constraints can be improved. We 
incorporate recent LHC limits on SUSY particles using SModelS [85, 86] - a tool that 
exploits the simplihed models results of the SUSY searches from the ATLAS and CMS 
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collaborations. Finally, after characterizing the remaining parameter space, we discuss 
the different channels available to further probe light stop scenarios in the MSSM. These 
include LHC SUSY searches, the associated production of stops with light Higgs, the 
monojet search for degenerate stop - neutralino, searches for heavy Higgses as well as 
direct/indirect DM detection. 

This paper is organized as follows. The set-up for the analysis is presented in 
Section 2 followed by a summary of the indirect and direct constraints in Section 3. Our 
results for the remaining available parameter space of the pMSSM with stops below 1.5 
TeV are discussed in Section 5 together with the potential for further probing the model 
with various collider and astro-particle observables. Our conclusions are presented in 
Section 6. 

2 Analysis set-up 

Here we have chosen a simplihed version of the pMSSM where only the ten parame¬ 
ters most relevant for the Higgs and DM sector are let to vary. All other parameters 
(squarks of the hrst and second generation and all sleptons) are hxed to a value large 
enough to evade all the LHC constraints. The ranges of the values of pMSSM param¬ 
eters we consider are listed in Table 1. We performed a flat random scan for values 
of parameters in these ranges. Note that sleptons close in mass to the LSP can give 
an important contribution to DM (co-)annihilation, for example this occurs frequently 
for staus within the constrained MSSM. We ignore this possibility here since the im¬ 
portance of coannihilation will be illustrated with third generation squarks, moreover 
coannihilation with sleptons rely on hne-tuning of parameters from unrelated sectors 
when no assumption is made about the underlying high scale model. 

We have used SuSpect-2.41 [87] to compute the pMSSM mass spectrum for a given 
set of input parameters. This includes two-loop corrections to the Higgs mass and NLO 
corrections to SUSY particle masses. In the hrst step of our scan, we have generated 
~ 0.75 Million pMSSM points for which the following statements apply, 

• All the criteria in SuSpect-2.41 for a theoretically valid point are satished, 

• The pMSSM spectrum has a neutralino LSP, 

• The lightest CP even Higgs boson mass (as computed by SuSpect-2.41) rrih sat- 

ishes llSGeV < < 130GeV, 

• The lightest stop and the LSP satisfy < 1500 GeV and m^o < 800 GeV 
respectively. 
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Parameter 

Scan range 

U(l) gaugino mass parameter: Mi 

20 - 2000 

SU(2) gaugino mass parameter: M 2 

100 - 2000 

Ratio of the vacuum expectation values of the two 
Higgs doublets: tan/3 

2-55 

Higgsino mass parameter: p 

100 - 3000 

Pseudo-scalar mass parameter: tua 

100 - 2000 

Stop tri-linear coupling: At 

-5000 - 5000 

Sbottom tri-linear coupling: Ab 

-5000 - 5000 

Mass parameter for the left handed third generation 
squark doublet: uIq^ 

100 - 2000 

Mass parameter for the right handed stop: rrijj^ 

100 - 2000 

Mass parameter for the right handed sbottom: 

100 - 2000 


Table 1. Scan ranges for the pMSSM parameters, all dimension full parameters are in GeV. 
The values of all the slepton mass parameters as well as ^^ > ^Di 2 
to 2 TeV. All the A-terms other than At and Ab are assumed to vanish. 


Once the above step is done, we assess the impact of Higgs boson signal strengths, 
flavonr violating observables, DM relic density and direct detection cross section, LEP 
data and finally, direct SUSY searches at the LHC on the selected set of ~ 0.75 Million 
pMSSM points. 

Note that we have checked the vacnnm stability and the absence of charge and 
colonr breaking minima via SnSpect. A more refined analysis along the lines of [88-90] 
might improve the constraints in the large At region. 

3 Indirect constraints 

• Higgs data : In order to stndy the compatibility of the pMSSM models with 
Higgs data we nse HiggsBonnds-4.1.0 [91, 92] and HiggsSignals-1.1.0 [93] which 
are linked to FeynHiggs for the compntation of Higgs mass and signal strengths. 
As far as the theoretical nncertainty in the Higgs mass calculation is concerned, 
we use the estimate given by FeynHiggs-2.10.0 [94]. While calculating the p-value 
in HiggsSignals, we set the number of free model parameters Np = 10. The p- 
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value is required to be more than 0.05 for an allowed parameter point. Note that 
the range for the Higgs mass assumed in our preselection will automatically be 
reduced by imposing these constraints. 

• Flavour data : Flavour physics has played a crucial role in the construction of 
the SM as well as constraining NP beyond the SM. The flavour structure of the 
SM is indeed very special, and any generic NP model suffers from large flavour 
violations in contradiction with the wealth of data from B-factories and recently 
also from LHCb. In the pMSSM, because of the degeneracy of the hrst two squark 
generations the flavour violation involving the hrst two generation of fermions is 
mild. Moreover, as we have set the masses of the hrst two generation of squarks 
to 2 TeV (which is allowed by direct searches at the LHC), their contributions 
decouple. Thus, the havour constraints mainly arise from processes involving 
the third generation of quarks, for example, decays involving 6 —)■ s transition 
In particular, the havour changing H-meson decays e.g., the radiative decay 
B{Bd —)■ Xs'j) and the fully leptonic decay B{Bs —)■ are known to place 

important constraints on the MSSM parameter space. In this work we have used 
the following limits for these branching ratios. 


2.78 X 10"* < X.i) < 4.08 x 10"*, 

1.43 X 10“’ < B(B, -> < 4.37 x 10“®. 


(3.1) 

(3.2) 


In addition, we have also imposed the following limits on the branching ratio of 
Bd —t , 


0.79 X 10"^° < B{Bd < 6.80 x 10"^° . 


(3.3) 


In the numerical analysis, we have used the public code SuperIso-3.3 [103] to 
compute these branching ratios. 

Note that we have not considered a few other potentially important observables 
such as the anomalous magnetic moment of the muon, — 2)^, the branching 
ratios of the two body leptonic decay of the H-meson B{B ^ ru) and the three 
body semileptonic decay B{B —)■ D{D*)tp). As far as {g — 2)^ is concerned, 
the SUSY contribution needs to be at least 1.0 x 10“® in order to be consistent 

^Interestingly, the LHCb collaboration has reported hints of NP in some of the B meson decay 
modes involving quark level b ^ s transitions [95, 96]. Although, there are pending issues with the 
reliability of the theoretical SM predictions and firm claims of the existence of NP can not be made 
yet, several NP explanations of these “deviations” have been proposed [97-101]. Unfortunately, a NP 
explanation within the MSSM seems unlikely [102]. 
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with the measured value at the 2cr level [104, 105]. This requires the existence 
of light sleptons and electroweak gauginos which can be easily achieved in the 
pMSSM [53]. We set the slepton mass parameters to a high value for simplicity. 

The leptonic decay B ^ tv which has tree level SUSY contribution from the 
charged Higgs exchange diagrams is also known to provide stringent constraints 
on the SUSY parameter space [106]. Interestingly, the most recent measurement 
by the Belle collaboration has brought down this branching ratio to a much 
smaller value than earlier measured [107], thus relaxing the tension with the SM 
prediction. All our points satisfy this constraint since many other constraints 
force the charged Higgs to be rather heavy any way. 


Let us now discuss about the three body semileptonic decay ^ Dtv) and 
B{B ^ D*tv). The BaBar collaboration measured the two quantities [108] 


R(D) = 


B{B —)■ Dt Vr 


and R{D*) = 


B{B D*t-iJ^) 


B{B^D£-v^) “ ^ B{B^D*i-ve)' 

and reported a 3.4a deviation from the SM when the two measurements are 
taken together. This result motivated a number a phenomenological studies 
both in the context of specihc models [109] as well as model independent ap¬ 
proaches [110, 111]. The BaBar collaboration itself ruled out a Type-H Two Higgs 
Doublet Model (THDM) at 99.8% conhdence level for any value of tan (3/mH± 
based on this data. The same would apply to the Higgs sector of the MSSM 
which is also a Type-H THDM at the tree level. However, the existence of two 
neutrinos in the hnal state of these decays makes their measurements quite chal¬ 
lenging experimentally and a conhrmation of these results by another independent 
experiment (e.g.. Belle H) is awaited. 


• Dark matter relic density and direct detection : The DM relic density has been 
measured precisely by PLANCK, = 0.1192 ± 0.00024 at 68% CL [112]. We 
impose the following upper bound 

Qh^ < 0.129 (3.4) 

which corresponds to the measured value after adding a 10% theoretical uncer¬ 
tainty - this number is a rough estimate of uncertainties that can arise for example 
from one-loop corrections to DM annihilation cross section [113-117] We im¬ 
pose only the upper bound to allow for the possibility that the neutralino is only 
a fraction of the DM. 

more precise value for the relic density has been released recently by PLANCK [118], this 
however has no impact on the results presented here since the theoretical uncertainty we assume is 
dominant. 






We also impose the mass dependent npper bonnd on the WIMP direct detection 
cross section obtained by the LUX experiment. For this, we have fitted the LUX 
upper bound [119] to an analytic formula which is given by, 


logic ^si^^ = 


7.029 


7.161 


8.569 


(log 10 WIMP logiomwiMP 
+ 0.755 log^o i^wiMP — 0.003 (log;^o i^wimp)^ 


Based on this we apply the following constraint, 

^ / C^LUX 

<7SI < ^CTsi 

where 


(3.5) 


(3.6) 




1 if 0.1103 < < 0.1289 , 


0.1196 


if < 0.1103 . 


(3.7) 


That is, if the relic density computed assuming the standard cosmological scenario falls 
below the PLANCK range, we consider that the neutralino constitutes only a fraction 
of the DM and explicitly ignore the possibility of regenerating DM although we will 
comment on this assumption in section 7. Moreover we do not make any assumptions 
about what would constitute the rest of the DM. Note that DM can also be searched 
for by indirect detection. In our analysis we have neither imposed the constraints 
from FermiLAT on photons [120], from PAMELA on antiprotons [121] as well as the 
preliminary limits from AMS on antiprotons [122], nor have we made any attempt 
to explain the anomalies observed by PAMELA [123] and AMS [124] on the positron 
spectrum. We do however, briefly discuss the impact of indirect searches in section 7. 
We use micrOMEGAs-3.5.5 [125] to calculate the DM relic density as well as the direct 
and indirect detection cross sections. 

In the MSSM, it is well known that the composition of the neutralino is crucial for 
determining the DM properties. For neutralino annihilation to be efficient enough to 
have VLh? < 0.129 requires either a LSP with a signihcant Higgsino or wino fraction 
or special tuning of parameters. In fact a dominantly Higgsino/wino LSP with a mass 
in the range 80 GeV to 1 — 2 TeV typically leads to < Dplanck^^ because 
of efficient annihilation into W pairs. A mixed state with some bino component is 
therefore preferred. Direct detection cross section on the other hand is large for a mixed 
gaugino/Higgsino LSP, in particular for the one that leads to the exact range of the 
relic density determined by PLANCK. The only possibility for such mixed neutralino 
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would be to lie above the TeV scale where the direct detection limits are weaker, we 
do not consider these masses since we want to highlight the SUSY spectrum below the 
TeV scale. Both cases with pure wino or pure Higgsino DM easily evade the direct 
detection constraints although the DM relic density cannot be entirely explained by 
neutralinos. Another possibility which allows also for a dominantly bino LSP consists 
in adjusting parameters such that ~ mz/‘^ or rn/i/2 or mu/2 thus providing a 
resonant enhancement of the cross section or having my ~ mg fermion- The contribution 
of coannihilation channels then reduces the relic density. The dominantly bino LSP is 
only allowed for light sfermions or when the mass is such that one can beneht from 
annihilation through a resonance in s-channel. Thus we expect to hnd a large number 
of scenarios with dominantly Higgsino or wino LSPs. 

4 Direct search constraints 

• LEP limits ; 

The generic limits from LEP are obtained directly from micrOMEGAs and mainly 
exclude charged particles. The lower limit on chargino is 103 GeV while those on 
sleptons, in particular staus are slightly weaker. In addition we have also imposed 
an upper limit on the Z invisible width, P^ < 2 MeV as well as constraints on 
neutralinos from a{e~^e~ —>■ < 0.1 pb where the heavy neutralino, 

decays mostly into hadrons and a LSP [126]. Such analyses constrain the very 
light neutralino LSP region, which is also strongly constrained by the upper bound 
on the DM relic density. 

Finally we impose the condition that all particles decay promptly (cr < 0.05 m). 
A strict requirement is that charged particles are not long-lived at the cosmologi¬ 
cal scale, we impose this more restrictive criteria because the LHG limits that we 
implement below assume prompt decays leading to Missing Transverse Energy 
(MET) in the hnal state. 

• LHG limits : Searches for SUSY at the LHG form an important ingredient to 
assess the viability of the scenario under consideration. Here, we describe our 
procedure for evaluating the LHG constraints. We use SModelS [85, 86], a tool 
designed to evaluate the LHG constraints on NP using simplihed model spectra 
(SMS) results. SModelS is designed to decompose the signal of any arbitrary NP 
spectrum with a Z 2 symmetry into simplihed model topologies and test it against 
the existing LHG bounds in the SMS context. The input to SModelS can either 
be an SLHA hie [127] containing the SUSY production cross sections a and the 
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branching ratios for the SUSY decays, B or a. LHE hie [128]. For this work we used 
the SLHA input containing a and B. The format for writing the production cross 
sections is specihed in [129]. The production cross sections are computed using 
Pythia6.4.27 [130] and NLL-fast-2.1 [131-138]. Given the information on the 
a and B, SModelS computes a x B for each possible decay of SUSY particles. 
The information of relevance to check the results against the LHC limits is the 
mass vector of the SUSY particles, the SM decay products and the a x B of 
the resulting topologies. A topology resulting from such SLHA decomposition is 
considered if the a x B > cxcut, with cxcut set to 0.01 fb. 

When dealing with an arbitrary NP spectrum, care must be taken to identify 
regions of compressed spectra as the decay products in such cases are not detected. 
SModelS ignores such soft decays when the mass gap between the mother and 
the daughter particles is less than the user dehned minimum mass gap, here we 
take 5 GeV. The resulting a x H for various SMS topologies were tested against 
ATLAS SMS interpretation for searches [139-151] and CMS searches [152-164]. 

5 Results 

We start our discussion with Table 2 where the impact of the various experimental 
data on the pMSSM parameter space is shown. Since we are randomly scanning over 
10 uncorrelated parameters, it is important to generate enough points to populate all 
dimensions. The table shows the effect of each successive experimental constraint on a 
well populated flat scan over 10 dimensional parameter space. The number of pMSSM 
points which survive after each successive constraint is presented in the second column. 
Although no statistical meaning can be attached to these numbers. Table 2 and the 
accompanying discussion helps to get an understanding of how the various type of 
data constrain the MSSM parameter space. It can be seen that out of the ~ three 
quarter of a million points that were generated, only 60% of them successfully satisfy 
the observed Higgs boson mass, signal strength data and other Higgs boson searches 
(which are implemented in HiggsBounds-4.1.0 and HiggsSignals-1.1.0). The severe 
impact of the measured branching ratio of Bd —)■ Xgj is also clear from the table. 
Out of all the models which satished Higgs data, only about 50% survive after the 
consistency with —)■ Wy branching ratio is imposed. It is worth mentioning that 

there is some amount of tension between the Higgs mass and the branching ratio of 
Bd — Xs'j in the MSSM. The latter gets contributions from chargino loops as well as 
charged Higgs loops. As the charged Higgs loop interferes constructively with the SM, 
its contribution is always positive. The chargino contribution on the other hand can be 
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Constraints 

No. of models 

Theory -|- neutral LSP -|- 

llSGeV < irih < ISOGeV -|- < ISOOGeV -|- m^o < SOOGeV 

741605 

Indirect bounds 

-h HiggsBounds -|- HiggsSignals 

435021 

+ —)■ Xg'^) 

211313 

+ B{Bs^d ) 

177961 

+ nhl + 

111167 

Direct bounds 

+ LEP 

59425 

-h Long-lived chargino 

30754 

+ LHG 

29266 


Table 2. Number of surviving pMSSM models after each cut. 


both constructive or destructive. For example, the stop—Higgsino contribution to the 
amplitude is proportional to mf nAt tan 13/m~ and hence, this contribution is enhanced 
for a large At. On the other hand, consistency with the measured value of the Higgs 
mass requires a large At if a light ti is desired. 

The DM constraints also reduce signihcantly the number of allowed points. In 
particular when the LSP is lighter than the W the only allowed points are near mz/2 
or m/i/2 corresponding to the annihilation through a Z or Higgs resonance. For heavier 
LSP masses the relic density upper bound basically selects LSP with large Higgsino or 
wino LSP, barring special conhgurations where coannihilations are important. Further¬ 
more the mixed bino(wino)/Higgsino can be in conflict with the LUX direct detection 
bound as mentioned above. We therefore expect that the LSP will most of the time be 
dominantly Higgsino or wino which in both cases implies that it is almost degenerate 
with the chargino. Such degenerate chargino - neutralino lead to long lived charginos, 
leading to charged tracks at the LHC. As explained before, these scenarios are not con¬ 
sidered in the present study and hence nearly half the points are discarded. Similarly 
the LEP limits on charged particles also rule out nearly half of the allowed points. 
Finally, LHC limits from the 8 TeV run rule out about 5% of the allowed points. Be¬ 
low we discuss in more details the impact of LHC constraints from SUSY searches as 
obtained with SModelS . 

In Fig. 1 we plot the allowed and excluded points after applying SModelS . In 
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m-,^ [GeV] ni'i^ [GeV] 


Figure 1. Summary of the allowed and excluded points by SModelS . In the left panel, the 
excluded points are displayed on top of the allowed points while in the right panel the allowed 
points are on top of the excluded points. 


the left panel we plot the excluded points on top of the allowed points while in the 
right panel, the plotting order is inverted. Clearly many points with light stop masses 
are not excluded by SMS results as implemented in SModelS -1.0.3. It is possible to 
exclude many points up to the maximum stop mass considered (1500GeV). However, 
no SMS result has reach for LSP masses greater than 300 GeV and therefore the region 
with higher LSP masses remains unconstrained. 

The right panel shows that different conhgurations of the MSSM spectra can evade 
the SMS results, thus allowing very light stop masses — even below 200 GeV. Despite 
the fact that SModelS combines topologies with the same hnal states and similar 
mass vectors, the main reason for the allowed region is that the SMS results obtained 
by LHC collaborations and used by SModelS assume a 100% branching ratio for the 
decay under consideration while in the MSSM branching ratios are often below 100%. 
For example consider a point in our scan with ~ 380 GeV and the rest of the spectra 
too heavy to contribute to limit setting. The dominant branching ratio for the stop is 
into (77%), however this channel cannot be exploited as the chargino decays into a 
virtual W and the LSP, a channel which is not implemented in SModelS. Therefore only 
the channel txi (with a HR 0.13) can be used to constrain the stop. This point gives 
a X B ^ 8.29 fb, whereas the experimental upper limit in CMS-SUS-13-011 derived 
assuming a 100% B is 438 fb. Hence this point does not get excluded by SModelS . 
The right panel of Fig. 1 thus justihes the possibility that a light stop is consistent 
with the current LHC SUSY searches as well as the observed Higgs properties and the 
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heavy Higgs searches, along with the DM direct detection limits. 
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Figure 2. Breakdown of the most constraining analysis for the points exclnded by SMod- 
elS in the — m^o mass plane (upper panel) and in the m~± — m^o mass plane (lower 
panel). For comparison, exclusion curves from various ATLAS and CMS stop searches are 
overlaid (solid lines). The dashed line shows exclusions arising from ATLAS searche for stop 
decays to cy? and to &Xi//h 


Figure 2 shows the breakdown of the most constraining analysis for the excluded 
points in the stop - LSP mass plane and chargino - LSP plane. For each of the excluded 
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points, we select the most constraining analysis It is possible that a point is excluded 
by more than one searches at the LHC. We overlay the exclusion lines obtained by 
ATLAS and CMS from stop searches to guide the eye. Indeed we see that most of 
the points excluded by the t —)■ txi searches fall within the corresponding exclusion 
contours while higher stop masses are in fact excluded by constraints coming from 
sbottom searches. There are also signihcant constraints arising from sbottom searches 
for the light stop masses in the regime where t —)• txi is kinematically forbidden while 
b — bxi is allowed. This kinematic conhguration is common when ti, bi are left-handed, 
thus sbottom searches provide indirect constraint on the light stop scenarios otherwise 
elusive at the LHC. For the kinematic edge where m^ — m^o < rrit, four body decays of 
stops as well as decays via cXi are also utilized at the LHC. Results from ATLAS [165] 
(for four body stop decay as well as decay via charm) and CMS [166] (stop decays via 
charm) searches are available. The ATLAS search yields stronger limits. Unfortunately, 
the cross section 95% C.L. observed upper limit map on the a x B is not available 
in [165]. This search hence could not be included in SMoDELSand was not used in 
this study. The red dashed exclusion lines obtained from the ATLAS searches [165] 
are overlaid in Fig. 2 (upper panel) for comparison. It can be seen that they do not 
cover a large region of parameter space, hence do not affect our conclusions drastically. 

Figure 2 (lower panel) shows the exclusions in the x^ — Xi plane, three distinct 
branches can be seen. The points along m^o 45 GeV or m^o ^ 60 GeV correspond to 
mostly bino-like y?, while the points along the diagonal represent dominantly Higgsino 
or wino Xi ■ The hgure illustrates that most of the times, the chargino is nearly mass 
degenerate with the LSP and either decays invisibly or via an off-shell W to LSP, thus 
evading the SMS limits which do not include either of these channels. The current 
searches for X 2 xf decaying to WZ -|- MET therefore exclude only a few points in this 
parameter space. These are located along the bino branches of the plot. Furthermore, 
along the diagonal lines where the chargino decay is invisible, the searches for stop and 
sbottom pair production with direct decays to LSP or Next-to-Lightest Supersymmetric 
Particle (NLSP) contribute the most to the exclusion. 

In Fig. 3 we plot the allowed points in the plane. We separate the total 

number of allowed points in three categories i) pMSSM points which have the lightest 
stop as the NSLP (red stars) ii) lightest sbottom as the NLSP (blue triangles) hi) 
lightest chargino as the NLSP (green circles). We have shown these three set of points 
separately to emphasize that most of the allowed parameter points where ti is the NLSP 
(the left panel in Fig. 3) lie close to the = m^o line. The requirements on the relic 

^The most constraining analysis is defined as the analysis which leads to the largest ratio of the 
theory cross section to the experimental upper limit. 
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Figure 3. Allowed points in the — m^o plane when the NLSP is or (left panel), 
m~± (right panel). The exclusions curves from the analyses of the stop searches used in testing 
the points in the — m^o plane. The lines with different colors and styles correspond to 
different SMS results for direct stop pair production decaying to tt + MET. 


density leading to the stop co-annihilation region is responsible for this strip. In these 
kind of scenarios, it is difficult to constrain these points via direct stop searches. The 
right panel shows that most of the points belong to the category of chargino NSLP 
as argued above. In such cases, the stop has cascade decays via chargino or heavier 
neutralinos, thus reducing the branching ratio for each mode and leading to weaker 
exclusions than expected from SMS results. Moreover only results where the chargino 
decays through a real W are included in SmodelS. 

We have demonstrated the exclusions obtained with the help of SMS results. How¬ 
ever, it might be possible to obtain stronger exclusions by means of recasting an anal¬ 
ysis, for example using the approach in [167-169]. This is clearly beyond the scope of 
this study. It is worth noting that allowing for light sleptons in the scan might lead 
to further exclusions driven by chargino - neutralino decays via intermediate sleptons. 
Furthermore, there are SMS interpretations for the decays of heavier stop e.g. stop 
searches with hnal states involving Higgs or Z boson [170]. As will be demonstrated in 
section 6.3 due to the constraint on the Higgs mass the t 2 is always heavier than 700 
GeV, where these searches currently do not have sensitivity. 

6 Probing light stop scenarios at the LHC 

From the results just presented it should be clear that the light stop scenarios in the 
pMSSM offer a variety of signatures at colliders. Note that this includes not only those 
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from direct stop production but also from other light super particles. In particular, a 
light left-handed stop means, quite often, a light left-handed sbottom as well. Thus, 
in this section we hrst investigate the main signatures that could not be constrained 
by SModelS and suggest additional topologies which may be pursued at LHC 14. Here 
we concentrate on hnal states produced by stop and sbottom decays and show missing 
topologies with large cross-section. We further also present missing topologies that 
arise from the electroweak sector for our allowed MSSM points. Then we examine 
other potential signatures from light stop associated production with an extra jet or a 
Higgs. Furthermore, we examine possible hnal states resulting from the decay of the 
heavier stop and we also discuss aspects of Heavy Higgs phenomenology for the set of 
MSSM points which are allowed in our light stop scenario. 

6.1 Improving simplified models interpretations at 8 TeV 

Figure 1 demonstrates the existence of a large number of points not excluded by SMod¬ 
elS . In order to understand the characteristics of these points and to suggest further 
ways to constrain the non-excluded regions of parameter space, it is interesting to ask 
which SMS topologies not covered by current SMS searches in the SModelS database 
prevail in these points. These are dubbed ‘missing topologies’ and will be discussed in 
this section in details. The missing topologies are derived purely on the basis of cross 
section times branching ratio {a x B) computations and hence do not take into account 
the sensitivity of the experimental searches. For example, in the results that follow, the 
hadronic decays of W will show up most of the times because the branching ratio for 
W decays to quarks is higher than into leptons, however, it is more difficult to beat the 
backgrounds while searching for hadronic decays. Note that it is possible to constrain 
some of the missing topologies by means of reinterpreting the existing experimental 
searches, thus missing topologies are not always associated with a new signature. 

The procedure used to derive missing topologies in SModelS is as follows. Each 
SUSY point leads to more than one missing topology. SModelS sums over the a x B 
for all the topologies irrespective of the mass vector of the SUSY particles. Up to 10 such 
topologies with the highest a x B are recorded in the output. In the following results, 
we suppress pair production of LSP which occurs as a missing topology. Moreover, we 
do not consider any initial/hnal state radiation effects. 

In Fig. 4, for each non-excluded point, the SMS topology leading to the highest 
a X B is chosen. We plot the seven most frequently occurring topologies originating 
from decays of stop/sbottom (upper panel) or electroweakinos (lower panel). The 
labels of the plot are written in a simplified notation with respect to the one used in 
SModelS [85, 86], however it is easy to map the current notation to the original one. 
In this work, every branch is enclosed in parenthesis () and every vertex is separated 
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Figure 4. Most dominant missing topologies in the stop - LSP mass plane. For each non- 
excluded point, the SMS topology leading to the highest a x B is chosen. The seven most 
frequently occurring topologies are plotted, the frequency of topologies indicated in the legend 
are ordered (top to bottom) with the exception of the label ‘other topologies’. The top panel 
shows the seven most frequently occurring topologies originating from stop/sbottom decays 
while the lower panel represents those originating from electroweakino decays. 


with a comma (,). In both notations, it is assumed that every branch is accompanied 
with MET at the end of the cascade. The possible origins of these topologies are 
explained in Table 3. Note that with this notation we do not distinguish particles 
and antiparticles and a sum over light quarks and light leptons is understood. In 
principle, decays apart from those illustrated in the table can contribute to the missing 
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topologies however, because SMoDELSdoes not keep track of the SUSY particles but 
only of their masses, this information is lost in the process of decomposition. Let us 
hrst concentrate on the upper panel of Fig. 4. The most frequently occurring topology 
is (b), (t), meaning b + t + MET. The topology occurs via asymmetric decays of 
pair produced stop (sbottom) with one of the stops (sbottoms) decaying directly to top 
(bottom) + LSP while another stop (sbottom) decays to bottom (top) via chargino. 
The decay of chargino itself is invisible when the chargino - LSP mass gap is less than 5 
GeV. Given the extreme degeneracy of the chargino - LSP masses the frequency of this 
topology is hardly a surprise. In general, this topology is difficult to be constrained by 
the current LHG searches as it leads to 2 b jets + one lepton + MET or 2 b jets + 2 jets 
+ MET final state. This is a topology with low jet multiplicity while most of the current 
searches for stops (sbottoms) require high multiplicity of light jets. Recently, limits for 
this mixed topology were made available by both the ATLAS [171] and GMS [172] 
collaborations. These limits assume exactly 50% branching ratios for each decay mode. 
The limits thus obtained are a sum over symmetric and asymmetric decay modes and 
do not represent upper limits on b + t + MET final state alone. For this reason, these 
limits are not applicable to our topology. 

Other topologies involve the cascade decay of one or both pair-produced squarks 
into a quark and a heavier neutralino (chargino) which then decays into a LSP and 
an off-shell Z (W). It is important to notice that the topologies (t) (t, , 

(t. WLd/Z*had^(t, (b, (t) for which stops contribute, are 

suppressed for > 1 TeV. This is simply because the production cross section of 
stops heavier than 1 TeV is extremely small. Thus the only relevant topologies are 
those resulting from sbottom decay, it occurs only when the right handed sbottom is 
lighter than the stops. Note that when or Z^^^ are found in a missing topology, 
there is also the same topology with Wi*^p or Z^^^ with a cross section reduced by the 
relative leptonic to hadronic branching ratio of the gauge boson. Despite the smaller 
cross section the leptonic final states typically have a much better signal to background 
ratio. 

The most frequent topology resulting from electroweakino production and decays 
represented in the lower panel of Fig. 4 corresponds to associated production of the LSP 
with a chargino (heavier neutralino) decaying via an off-shell W (Z). Other topologies 
get their dominant contributions from the decays of at least one heavier neutralino or 
chargino, as detailed in Table 3. The very mixed nature of these electroweakino lead 
to a sizeable production cross sections despite heavier masses. 
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Topology 

Decay 

(b)(b. WUZ’hJ 

bill bxi bxl^ bx\ b Z^^^Xi 

(b)(t) 

ti ti txi b Wsoft Xi 
bill bx\tWsoftx\ 

(b; bF^ad/-^Ld)(b) ^had/^had) 

till ^b Xi b bbjCad Xi 


till txi txl tx\ tZl^^Xi 

hli txi txl tzi^dxl tzi^dxl 

(b)(t. WUZkJ 

hli txl txl txl tZhadxl 

hli txl txl tzi^xl tzi^^xl 

(b ^w/-^Ld)(b ^had/^had) 

h h txl txl tzi^^xl tzi^^xl 
hh txl txl tW^adXltWi^adXl 

(b, WUZUit) 

Ml ^ bxl txl bZl^^xl tWsoftxl 
hli -)■ bxl txl bWhadXl txl 

(b)(b, z;^) 

Ml ^ bxl bxl bxl bZ^pXl 

(b)(b, 7) 

li li -t bxl bxl-^ bxl b XXI 
lih bxlbxl bZl^xlblxl 

(inv) 

xlxl W^adXlxl 
xlxl zi^dxlxl 

(W)(W) 

xlxl -^WxlWxl ^W.xlWWsoftxl 

(inv) (photon) 

xl X 2 ^ Wsoft xl 7 xl 

(inv)(W) 

xl ^ Wsoft xlWxl^ Wsoft xl W Wsoft xl 

(invXW'-,:,) 

xfx“,^ IV,Infix’, 

(inv)(^rep) 

xt x'i -X W'W. X? ZlX'i 


x'i^ ^ zi^^lw;,Xl 
xfxf ^ fVL^xl fv,:,x1 

(inv)(b, b) 

%xt ^ z;,%iv,,„fj 


Table 3. Missing topologies represented in Figs. 4,5 and 6, written in SModelS notation and 
the corresponding physical process. Wgoft represents soft decays of W, which are undetected. 

and hF;*p represent the sum over all three generations of leptons, and represent 
the sum over hrst two generation quarks. 
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The main topologies, mono - W ((inv)(W)), mono - photon ((inv)( 7 )) and 
diboson ((W)(W)) ^ cover different mass range s for the LSP and this independently 
of the mass of the stop. Mono - W topologies could be used to probe the region 
100 < m^o < 200 GeV while the mono-photon dominates for higher LSP masses. Note 
that the mono - photon topology is here associated with xf + X 2 production where 
the chargino decays invisibly and the heavier neutralino decays via a loop-induced 
decay into the LSP and a photon. Because of the small difference between the X 2 ~ Xi 
masses, the loop-induced decay can have a large enough branching ratio (typically 
0{5%)) to give a signature while the products of the three-body decay are too soft 
to be detected and are thus registered as a pure missing energy signature. Note that 
these photons could have a large pt- The diboson topology which arises from heavier 
chargino decays cover the full mass range, in particular the bino-like neutralino branches 
(with m^o ^ 60 GeV) where there are no topologies occurring from the decays of 
stops/sbottoms. Updates of searches like [143] in the diboson final state from Run-II 
will thus be important. In fact, reinterpreting some of the existing mono - lepton [173, 
174] and mono - photon [175] searches could be useful in further constraining the mono 
- W and mono - photon final states. This point is left for further investigation. Such 
channels could thus indirectly constrain scenarios with light stop masses beyond the 
current reaches of direct stop searches. Note that for all topologies involving virtual 
W’s or Z’s, it is the hadronic decays of the gauge bosons that dominate. However, the 
leptonic decays will also be present with smaller a x B but can in principle have higher 
sensitivity than the hadronic channels. 

As explained at the beginning of this section, each SUSY point leads to more 
than one missing SMS topology, a list of up to 10 such topologies is available from 
SModelS . In Fig. 4, the topologies with the highest cross sections are described. The 
question of which other topologies occur in this scenario, and whether they lead to 
some more sensitive final states still needs to be answered. In order to illustrate this, 
all topologies with a. a x B > Ifb are sorted according to their frequency of occurrence 
and the seven most frequent topologies are shown in Figs. 5 and Fig. 6. Once again, in 
order to make it easy to understand the origin and correlation of the topologies, only 
topologies originating from stop/sbottom decays are plotted in Fig. 5 in the — m^o 
plane, while those from the decays of electroweakinos are plotted in Fig. 6. Gomparing 
Fig. 5 to the top panel of Fig. 4 it is clear that along with the dominant topologies, 
(b) (t) and those involving hadronic W*/Z* decays that are found in both figures, it 

®In principle, SMS results for this topology are available by ATLAS [143]. However, this result is 
not yet included in SModelS . The reach of this analysis at the moment is very limited, excluding 
chargino mass up to 180 GeV and reaching up to maximum neutralino mass of 30 GeV. Inclusion of 
this search does not exclude any points as the LSP is always heavier than ^ 40 GeV for our scenario. 
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is also possible to find other topologies which might have a better signal to background 
ratio e.g. (b) (b, or (b) (b, 7 ). The latter again occurs from the loop-induced 

decay of the second neutralino produced in sbottom decay as explained in table 3. This 
topology with 2 b jets, 1 photon -|- MET in the final state can be used to constrain the 
kinematic edge with — m^o < 50 GeV. 

Missing topologies in Fig. 5 show a strong correlation with the sbottom mass. 
Topologies involving two b jets, MET and a virtual gauge boson originating from the 
decay of sbottoms into a heavier neutralino/chargino occur mostly in the region with 
a small difference between the sbottom and the LSP. Only a few topologies involving 
stop pair production extend beyond > 800 GeV where the sbottom pair production 
becomes too small. 
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Figure 5. The figure shows missing topologies with cross section greater than 1 fb originating 
from stop or sbottom decays and sorted on their frequency of occurrence. The seven most 
frequently occurring topologies are illustrated. Topologies originating from electroweakino 
decays are suppressed. 


In Fig. 5, it is difficult to highlight topologies originating from the decays of stops. 
The reason is the large mixing angle in the stop sector due to the Higgs mass constraints. 
In this case, the stop has no preferred decay channel, thus, missing topologies with a 
large a x B which naturally result from a large branching ratio in a single channel are 
absent. 

As the m~± and m^o are extremely mass degenerate, we plot the missing topologies 
in the electroweak sector in the X 2 ~ Xi plane. Fig. 6 shows that most of the topologies 
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Figure 6. The figure shows missing topologies with cross section greater than 1 fb originating 
due to electroweakino decays sorted on their frequency of occurrence. The minimum cross 
section of any topology plotted is Ifb. The seven most frequently occurring topologies are 
illustrated. Topologies originating from stop/sbottom decays are suppressed. 

in this sector are in fact off-shell decays of W and Z bosons associated with MET. 
Clearly, in most of the region where hadronic decays of W or Z are dominant, the 
leptonic decays are also present. Moreover the importance of the diboson missing 
topology is once again evident. Expectedly, in the region with > 400GeV, 

the leptonic decays become irrelevant due to low a x B. 

6.2 Stop NLSP at LHC 14 TeV 

In this section we focus on the region of parameter space where the lightest top squark is 
the NLSP and has a small mass difference with the LSP hence the decay ti —)■ txi is not 
open. If the mas gap mj_^ — m^o < 85 GeV (which is true for most of the allowed points 
with a stop NLSP) then the 3-body decay G —)■ 6 W~^ Xi is kinematically forbidden 
and the only allowed decay modes are the 4-body decay ti —t 6 x? / /' and the flavour 
violating decay ti —)■ c/uxi- Because of the small mass gap between the light stop and 
the LSP, the decay products are very soft on average making this scenario extremely 
challenging for the LHC searches. This explains the very low lower bound on the light 
stop mass > 275 GeV) obtained with the LHC 8 TeV data for this specific scenario. 
Hence, it is interesting to explore the prospect of the 14 TeV LHC for this scenario. 
It was shown that using monojet -|- large missing transverse energy final state a stop 
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mass > 300 GeV will not be ruled out even with 100 fb“^ integrated luminosity at 14 
TeV LHC [75]. The same question was also investigated in Ref. [80] using ax and Mt 2 - 
Assuming B{ti —)■ cXi) = 1 it was shown that the exclusion limit can be extended to a 
maximum of ~ 450 GeV (depending on the mass gap between ti and Xi) with 100 fb“^ 
data. Moreover, the use of charm-tagging could prove very useful in the future [80]. 




Figure 7. Leading order cross sections for titij and titih at the 14 TeV LHC. Mad- 
Graph5_aMC@NLO version 2.2.2 [176] was used for the computation. The triangles represent 
points that are potentially excluded from ATLAS searches for ti —)■ c -|- y? and stop four body 
decays. 


As there are not many handles to suppress the SM backgrounds in the stop NLSP 
scenario, the radiation of a light Higgs from one of the stops in stop pair production 
could provide an additional handle to discriminate the signal. It was shown in [177, 178] 
that the associated production of a stop pair with a light Higgs could be large in some 
region of parameter space, basically due to a large titih coupling. However, when 
taking into consideration the mass of the Higgs which requires a contribution from the 
stops and maximal mixing, the titih vertex is suppressed, hence the titih cross section 
is expected to be quite low. Moreover this process depends on the parameters of the 
model whereas the cross section for titij production is expected to be much larger 
and depends only on In Fig. 7 we show titih (left panel) and titij (right panel) 
production cross sections at the 14 TeV LHG for all the allowed points with the light 
stop as the NLSP. The triangles represent potentially excluded points due to constraints 
coming from ATLAS H t c -|- and stop four body decays. The production cross 
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section for titih final state is at least an order of magnitude smaller than that for titij. 
This makes the titih channel much less promising than the titij hnal state. Although 
there can be considerable gain in the background reduction if the Higgs in the hnal 
state is tagged but the existence of irreducible background like pp —)■ hZ{-^ vv) (with 
a 14 TeV cross section ~ 150 fb) makes it extremely challenging. For dehniteness, 
in Table-4 we provide a few example benchmarks with different ti decay modes. We 
avoid a choice of benchmark point within the excluded regions of the existing ATLAS 
ti —)■ c X? stop four body decays. 




m~± 

I3{ti cXi) 

B{h ^ 6//'^) 

B{ti ^ bx^) 

a{titih){fh) 

305.3 

295.3 

820.7 

1 

0 

0 

7.3 

372.4 

364.0 

1097.0 

1 

0 

0 

3.7 

328.8 

301.2 

1033.6 

0.41 

0.59 

0 

2.5 

314.3 

305.1 

309.5 

0 

1 

0 

6.2 

308.3 

260.1 

264.0 

0 

0 

1 

12.0 

353.7 

319.8 

322.8 

0 

0 

1 

4.1 


Table 4. A few benchmark points to show specific examples of different stop decay modes 
and also the 14 TeV titih production cross section in fb. While for the first three benchmarks 
the lightest stop is the NLSP, for the final three points the lightest chargino is the NLSP. 


In the first two benchmarks ti exclusively decays to cXi and in the following two 
points ti exclusively decays to the 4-body final state. In the final two benchmarks ti is 
not the NLSP - the Xi and X 2 beneath ti in the spectrum. As a consequence, the 
ti exclusively decays to bxf- The chargino in this case will eventually decay to ff'xi 
giving rise to again a 4-body decay of ti but because of the existence of an on-shell 
chargino in the decay chain, the kinematics will be different. In summary, the stop pair 
production in association with an extra hard jet will be the most promising channel for 
probing the stop NLSP region with small stop-neutralino mass gap. The prospect of 
stop pair production process associated with a Higgs boson does not look encouraging 
mainly because of the small production cross section. We find that with the Higgs mass 
constraint, the stop composition is such that the titih coupling is small. This means 
that the titih cross-section is suppressed even for light stops, the suppression being 
dynamic and not so much due to kinematics. The existence of a Higgs in the final state 
may provide an additional handle to combat backgrounds, hence a more focused study 
may be worthwhile. As far as the decay of stop is concerned, there are three distinct 
categories where the final state objects and/or the kinematics are different. Hence, 
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dedicated searches at the 14 TeV LHC for each of them should be carried out. 


6.3 Decays of heavier stop 

Searches for the heavier stop could provide an alternative for probing the light stop 
scenarios, the main issue for exploiting these searches at the LHC remains the mass 
scale of the heavier stop. Due to the large radiative corrections required to achieve 
the correct Higgs mass in the MSSM either large At or heavy stop masses is required. 
A large value of At introduces a large splitting in stop masses and pushes the heavier 
stop above the TeV scale. In the left panel of Fig. 8 we plot the allowed points in 
the At — rrit^ plane, while in the right panel the same set of points are plotted against 
the mass difference between two stops Recall that At is related to the 

naturalness of the SUSY spectrum (see Eq. 1.4). Fig. 8 shows that the value of At is 
constrained to be approximately above 1.8 TeV. This, in turn, gives a lower bound on 
the hne-tuning parameter A > 50 (assuming Amess = 20 TeV) which amounts to ~ 2% 
tuning. 
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Figure 8. The allowed pMSSM points in the At — (left panel) and At— (m^^ - tuj^) 
(right panel) planes. 


An optimistic scenario with A ~ 50 leads to ~ 1 TeV and {frit^ — rritj ~ 500 
GeV (the region enclosed by the thick red line). Thus, an interesting outcome of our 
analysis is the possibility of a SUSY spectrum with a heavier stop mass around a TeV 
with a large mass gap (~ 500 GeV) with the lighter stop. The discovery potential of 
this scenario at the 14 TeV LHG will in general depend on the decay channel of t 2 
which in turn depends on the masses of the other SUSY particles. We will consider 
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two specific decay modes t 2 ^ tiZ and t 2 ^ tih which are particularly interesting and 
have in general large branching ratios (as we will show below). The couplings t 2 -ti-Z 
(in the limit when the mass splitting between t 2 and ti is large and Ma-,Mh ^ Mz) 
and t 2 -ti-h can be written as [177, 178], 
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The mixing angle between the left and right handed top squarks, 9^ and the mixing 
parameter are dehned through the mass matrices [177], 


M| = 


m 


LL 


+ m? mtXt 


X, = A- 


mtXt 

h 


m 


RR 


+ mf 


rriLL = 


Qs 


tan jS 
1 2 

+ (- — - sin^ 9w) cos 2/5 M| 


m^R = + - sin^ 9w cos 2(5 M| 


(6.3) 


2 2 I 

mj- = m. 5 — 
t 2 


^1,2 

sin 26 't = 


2mtXt 


m~ 


m~ 


^II + min T \l + Ar^Xf 

— , cos 26 *t = 


^LL - m%R 


t2 








t2 


Eq. 6.1 shows that the t 2 -ti-Z coupling will be large if Xt is large (which is required 
for the Higgs mass unless both stops are extremely heavy). The rrit enhancement can 
also be understood by noting that in the large energy limit (large mass splitting between 
t 2 and ti ) the Z boson can be replaced by one of the would-be goldstone boson using the 
goldstone boson equivalence theorem. The t 2 -t\-h coupling has slightly more structure 
and has both the F-term and D-term contributions. However, this coupling can also 
be quite large in some part of the parameter space. Hence, these two decay modes are 
well motivated from the fact that large X^ is required for the Higgs mass. Of course, 
the presence of other particles below the t 2 mass makes the picture more complicated 
and many other decay modes can contribute. In Fig. 9 we show the branching ratios 
into gauge and Higgs boson for all the allowed points. One can see that for a large 
fraction of points the sum of B{t 2 —t tiZ) and B{t 2 —t ti h) is quite large supporting 
our analytic expectation. 
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Figure 9. Left panel: branching ratio of t2 —?■ tih vs. that of t2 —>■ tiZ. Right panel: 
branching ratio of t2 —?• tiZ vs. that of t2 —>• tiW. The region only below the blue dashed 
line is physical. 


Hence, the pair production of t 2 with the subsequent decays t 2 ^ ti Z and ^2 —t H h 
can be a very interesting channel to look at. The 14 TeV prospect of this scenario was 
studied in [79] by one of the authors. It was shown that a 4-5cr signal can be observed 
for ~ 1 TeV t 2 with 100 fb“^ integrated luminosity. We refer the readers to [79] for 
further details. 

6.4 Decays of heavy Higgs 

In this section we turn our attention to the neutral heavy Higgs bosons H and A. 
An important result of our analysis is that the mass of the heavy neutral scalar H is 
constrained to be larger than ~ 450 GeV. Note that this lower bound arises only after 
imposing constraints from the light Higgs signal strengths and Heavy Higgs searches 
implemented in HiggsSignals and HiggsBounds . Moreover the A is quite degenerate 
with the H, the mass difference being almost always less than 5 GeV, which is less than 
one percent of the common mass. For such heavy Higgses, the widths are ~ O.OSmj^, ruA 
from SM fermions alone [179]. These widths can only increase when decay channels 
into sparticles are included. Thus the mass difference between H and A, is always 
comparable or smaller than their widths. It is well known that the agreement of the 
observed signal strengths with those expected in the SM, actually forces the global hts 

TT 

to the alignment region where |(/5 — ~ ■ As a result the gauge boson couplings 
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are severely suppressed for the H. Hence, even with the large mass of H and the 
enhancement factor in the VV decay width due to the decays in the longitudinal V 
bosons, the branching ratio into vector bosons is not above one percent. Recall that the 
AVV vertex is absent at tree level. We have also checked that the decays H ^ hh and 
A ^ Zh have branching ratios smaller than one percent for all our points. Thus, the 
only relevant tree level decay modes can be into standard model fermions and sparticles: 
the sfermions and the electroweakinos. Moreover the branching ratios should be similar 
for the H and the A. 

The large mass of the H means that the decay to the tt hnal state is now kinemat¬ 
ically allowed. As a result it can be the dominant decay mode for small values of tan f3 
and the B{H — tf) can be as high as > 80% for tan/3 ~ 5. With increasing values 
of tan/3, B{H —)■ tf) drops gradually and the H ^ bb decay mode starts dominating. 
These two branching ratios sum up to about 80% in large part of the parameter space. 
The remaining 20% is mostly taken up by B{H —rf) decays. The SM decay modes 
are however suppressed when decay channels into SUSY particles become important as 
will be discussed below. In the left panel of Fig. 10 we show the 14 TeV cross sections 
for H production (cxf^rpgy) times the branching ratio to the tt and tt hnal states for 
the allowed points for which crf^rpgy x B{H //) is greater than 1 fb. The production 
cross section were computed with SusHi-1.5.0 [180-188]. The decays of H/A into r pair 
clearly offer a search channel with crf^-pgy x B that can reach 0(10) fb, for tan/3 ~ 20 
and niH ~ 1 TeV. For the tt hnal state crf^rp^y x B{H ^ tt) can be as high as 10’s of 
fb and higher for low tan /3 ~ 15 and tuh ~ 1000 GeV. 

The decay of the heavy Higgs H to SUSY hnal states can also be quite important. 
In Fig. 11 we plot the product of cross section times branching ratio for H to the 
electroweakino hnal states All channels where the heavy Higgs can decay to charginos 
or neutralinos are summed over. The ‘invisible’ decays of the H include XiXi s-s well as 
the decay into heavier states when the mass diherence between the sparticle (eg. X 2 ) 
and the LSP is small; ~ 5 GeV. We hnd that the total ‘invisible’ branching ratio of 
the H is always less than 30% and mostly less than 10% while the visible decays into 
electroweakinos can reach 80%. The branching ratios for the A are also similar to those 
for the H presented in the hgure. In fact, apart from the production cross section, the 
phenomenology of H and A for these allowed points are rather similar. It is also well 
known that for the same mass, the gluon fusion production cross section can be higher 
for A than for H, the exact values depending on the masses of the squarks. 

There have been dedicated discussions of the Higgs sector of the MSSM, allowed 

®There is a small fraction of H that decays to light stops as well. However, this channel is not 
always kinematically open and most of the times B{H —>■ tit\) < 0.1. 
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Figure 10. The product of 14 TeV cross section of the CP-even heavy Higgs and its branching 
ratio to rr final state (left panel), tt final state (right panel) in the tan/? — Ma plane. 
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Figure 11. The product of 14 TeV cross section of the CP-even heavy Higgs and its branching 

ratio to electroweakino final state in the mn -tan /? plane (left panel) and in the M 2 - [i 

plane (right panel). 


after the LHC Run-I, beginning from the analysis of the partial hrst data set [189] 
in the pMSSM framework to the more recent analyses [179, 190] in the liMSSM and 
pMSSM scenarios respectively which focus on searches for the Heavy Higgses at the 
LHC 14 TeV. 
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Obviously, as the above discussion shows the tt hnal state is perhaps the most cru¬ 
cial in the low tan range and the rr in the large tan /3 range. While the tt channel 
has received a lot of attention in the past, the LHC Run-I results [191] have forced 
attention to be focused on the tt hnal state. It has been known since a long time that 
for this hnal state, interference with the tt QCD background gives a very characteristic 
peak-dip structure [192, 193]. The feasibility of using it to isolate the signal from the 
background as well as the diherence in the spin spin correlations between the t and i for 
the background and the resonant signal have been discussed in the literature [194-196]. 
In the CP conserving case, the H and A amplitudes do not interfere but still the pres¬ 
ence of almost degenerate H, A can degrade the ehect. The peak-dip structure seems 
to be subdominant to the ehect of higher order corrections [197] and more intricate 
cuts may have to be devised to enhance the resonant Higgs contribution. Moreover 
a recent analysis [198] concluded that this peak-dip structure will get degraded due 
to the limited resolution of the tt invariant mass, the statement being even more true 
with the presence of degenerate H and A as would be the case here. Hence kinematic 
cuts which exploit the ehect of the spin-spin correlations for the tt produced in the H 
decay would be necessary (see for example [199]). A simple analysis of [179] which 
includes such cuts, shows that at the LHC 14TeV one could be sensitive to tan ~ 6 
for ruH ~ 500 GeV and to tan /d ~ 1 for mn ~ 1 TeV. More analyses, to improve the 
sensitivity of this channel are required, see for example [198, 200]. 

Of some interest are the decays into the electroweakinos where the product of cross 
section times branching ratios can reach values as high as 10^ to 10^ fb for H masses 
up to 700 — 1000 GeV and tan /3 ~ 20 — 30. Thus we hnd that this channel can oher 
interesting search possibility for the H/A. Note that the general conclusion of [179] that 
the parameter ranges which gives large branching into these channels have been ruled 
out by the LHC trilepton constraints does not apply in our analysis where Mi, M 2 values 
are not related. The large values of the cross section times branching are concentrated 
in the low M 2 -/i region and hence the chargino/neutralinos produced are likely to give 
rise to hnal states with real or virtual WjZ. However, the topology for hnal states 
resulting from the H decay into electroweakinos, followed by their cascade decays, will 
depend on their couplings and on the masses which are strongly inhuenced by the relic 
density upper limit. A detailed investigation of the potential to use electroweakino 
decay channels to extend the reach for heavy Higgs, a topic that was hardly addressed 
[201] is beyond the scope of this paper. 

As mentioned already the H has an invisible branching ratio which can be upto 0.3 
but is mostly less than 0.1. One can in fact search ‘directly’ for such an invisible H at 
the LHC via the associated production of the H with a vector boson [202] or production 
of H via the vector boson fusion [203] . Indeed currently bounds exist on the invisible 
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branching ratio of the 126 GeV Higgs by both the CMS [204] and ATLAS [205, 206], 
using both the modes. For a 120 GeV Higgs boson, the 14 TeV LHC should be able to 
probe a value of the branching ratio as low as 0.17 [207]. Earlier projections [202, 203] 
had looked at larger values of the scalar mass. Unfortunately, these channels will be 
of not much use in the present case for the H as its couplings to a VV pair are highly 
suppressed. Hence the production of H in gluon fusion with associated jets followed by 
the decay of H in invisible channel [208] or associated production of H with a ti pair 
followed by an invisible decay of the H, are the two possibilities for such an invisibly 
decaying H. The estimates of the expected rates for the latter channel presented in [198] 
assuming an invisible branching ratio of 1, shows that the search in the channel ti pair 
+ MET will be a challenging one. As far as an invisibly decaying A is concerned, 
for the same mass and coupling strengths, the associated production will have smaller 
rates for ttA compared to tiH (see for example [209, 210]). Hence production of the 
(pseudo) scalar in gluon fusion with associated jets will be a better channel to probe 
for such an invisibly decaying A. Glearly, more studies are required. 

7 Complementarity with dark matter searches 
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Figure 12. Spin-Independent y^-nucleon cross section vs the mass for ti (left), bi (cen- 
ter)and xf (right) NLSP. For comparison, the expected limits from XENONIT and the 
neutrino coherent scattering are superimposed. 


The spin-independent (SI) neutralino-proton cross section is displayed in Fig. 12 
for the cases of the ti, bi and Xi NLSPs. The points are color coded according to the 
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stop-LSP mass difference. Although this quantity is not directly relevant for direct 
detection it is useful to highlight the complementarity with collider searches. Here the 
cross section is rescaled according to Eq. 3.7 to take into account scenarios where the 
neutralino is only a part of the DM. 

The results can be understood knowing that the predictions for the SI cross section 
are basically governed by the nature of the neutralino, pure states leading to small cross 
sections and mixed Higgsino/gaugino to the largest. A large fraction of the points will 
be probed by XenonlT, in particular those with a chargino NLSP - which constitute 
the largest sample - since they are typically associated with a dominantly, yet mixed, 
Higgsino(wino) LSP. The right panel of Fig. 12 shows that XenonlT has the potential 
to cover the vast majority of points where 40 GeV < — m^o < 85 GeV. Those 

are the points that are not well constrained by current LHG bounds because the stop 
decays mostly into hxi and the chargino in turn decays via a virtual W. Some of the 
points with — m^o < 40 GeV can also be probed by XenonlT provided the chargino 
is the NLSP. Furthermore the scenarios with a LSP with a mass around 45 or 60 GeV 
(corresponding to the so-called bino branches mentioned in previous sections) should 
be entirely probed with Xenon-IT. A few points with chargino NLSP he below the 
coherent neutrino scattering background, these are typically associated with a pure 
Higgsino or wino LSP. 

Less promising for direct detection are scenarios with a squark NLSP. Figure 12, 
left and center panels, shows that only a small fraction of the points will be probed by 
the future XenonlT detector, a few points even lying below the irreducible coherent 
neutrino scattering background. In particular some of the points where — m^o < 
40 GeV, in the left panel of Fig. 12, can lead to a very small cross section. The 
reason is that for such mass splitting the value of the relic density is governed by the 
coannihilation channels with stops, hence a dominantly bino LSP is allowed. Its weak 
coupling to the Higgs lead to a suppressed SI cross section. Thus these points that are 
hard to probe at the LHG can also evade direct DM searches. Note that going beyond 
the standard cosmological scenario by assuming that DM can be regenerated , for 
example from decay of moduli helds re-injecting neutralinos after the freeze-out [211- 
213], such that .^ = 1 would lead to much stronger constraint [214]. Many of the 
scenarios would already be constrained by LUX and nearly all of the ones with chargino 
NLSP could be probed at XenonlT. 

We have also computed the indirect detection cross section for LSP annihilation 
into T~^T~, bb and W~^W~ and compared this with the exclusion obtained by FermiLAT 
from observations of the photon flux from dwarf spheroidal galaxies (dSphs) of the Milky 
Way [120]. Fig. 13 shows the results for different channels after rescaling by Since 
a large fraction of the points have ^ < 1, the rescaled cross section is often strongly 
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Figure 13. Indirect detection cross section as a function of LSP mass in the WW, rr and 
bb final states . For comparison the current FermiLAT [120] limit is superimposed. 


suppressed, hence only a handful of points in the bb channel are excluded by such 
searches - basically those where the cross section is enhanced by annihilation through 
a heavy Higgs. Again assuming ^ = 1 would lead to a completely different picture, 
with a strong increase in the predictions of the cross sections. In particular av for DM 
annihilation in the WW channel exceeds the FermiLAT limit for most of the scenarios 
with a LSP mass below SOOGeV, in agreement with the results in [214, 215]. It was 
also shown that PAMELA limits from antiprotons can constrain such scenarios [216] 
and that the wino can be also effectively probed by FermiLAT searches for gamma-ray 
lines from neutralino annihilation into photon pairs (or jZ) [217]. 

8 Conclusions 

The discovery of a SM-like Higgs boson has strong implications for SUSY since a 
Higgs mass of 126 GeV requires large radiative corrections from the stop sector. This 
requires heavy stops and/or large mixing which is in conflict with the naturalness 
arguments. Stops therefore play a central role in SUSY and are a key ingredient in 
testing the naturalness of the MSSM. The LHG has performed many dedicated searches 
for stops. However, their production cross sections are small compared to the hrst two 
generation squarks, which leads to degraded LHG runl limits on their masses. Hence, it 
is interesting to investigate to which extent light stops are still allowed and demonstrate 
various possible future probes of the resultant MSSM scenario at the LHG via not just 
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the searches for the stops but also for other sparticles such as sbottoms, electroweakinos 
and even the heavy Higgs. 

We hrst determined the regions of the pMSSM with ten free parameters compatible 
with a light stop (specihcally with a mass below 1.5 TeV) after taking into account 
current constraints from the Higgs mass, Higgs signal strengths, flavour physics, the 
upper limit on the neutralino relic density as well as DM direct detection. Flavour 
constraints are very restrictive since they are in tension with the Higgs mass. For 
example, constraints from —?• can become stricter with large At, which is 

what the Higgs mass requires. The PLANCK upper limit on the relic density and the 
direct search limits from LUX are also in tension, combined they favour either TeV 
scale LSP or a almost pure Higgsino or wino LSP. The latter implies a supersymmetric 
spectra with a chargino NLSP and small mass differences between the chargino and 
the neutralino LSP. Another possibility is a squark NLSP (in particular a stop) since 
coannihilation can be used to obtain the measured value for the relic density. Both 
cases entail that LHC searches for SUSY are difficult. 

We have then used SModelS for all the MSSM points allowed by all the previously 
mentioned constraints to hnd implications of the limits obtained by the LHC collabora¬ 
tions on supersymmetric particles with the simplihed model framework. We have found 
that stops below 500 GeV can be consistent with all the LHC searches (including direct 
stop and sbottom production). The reason is either a compressed spectra or reduced 
branching ratios into the channels excluded with simplihed models. This is our hrst 
main result. 

Our results generally agree with those of [45] which also thoroughly investigates the 
impact of LHC SUSY searches on the pMSSM parameter space and found that light 
stops were allowed. Our approach however dihers in two ways. First in our choice of 
the set of free parameters. Most importantly, in contrast to us [45] allows light sleptons 
and gluinos. Second, we impose only an upper bound on the DM relic density leading 
to a large fraction of points with large higgsino or wino components. 

Another important aspect of our analysis is the identihcation of hnal states which 
would, in principle, be capable of probing the points surviving after the application of 
LHC, havor, dark matter and Higgs constraints. This analysis can help determine how 
best to extend the search reach at the LHC, particularly for the low masses still allowed 
after LHC constraints. One of the most important missing topology corresponds to 
asymmetric stop decays - for example one of the pair produced stop decaying to txi 
and the other into bxf with the chargino being invisible because it decays into the 
LSP and soft jets. Other possible signatures which might improve the reach involve 
the decay of squarks into a heavier neutralino or chargino which then decay to a LSP 
via virtual W/Z. 
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An extension of the reach in the region of small mass difference of the stop with 
the LSP entails using the jet and MET signature from a stop decaying into charm 
neutralino or via 4-body decays as done by the LHC collaborations. Although we have 
not used these channels as they played a minor role in our analysis, they should lead 
to strong constraints in the future. We also investigated whether the region near the 
kinematic boundary could be probed by considering associated production of stops 
with a Higgs. In principle tagging the Higgs allows to handle the background very well. 
We hnd that with the Higgs mass constraint, the stop composition is such that the 
titih coupling is small. As a result this cross section is small, even for small masses 
of the stops being considered, the suppression being dynamic and not so much due to 
kinematics. Hence the stop search in this channel is challenging and needs more detailed 
studies. The associated stop stop jet process is perhaps a better option, for the region 
near the kinematic boundary, though potentially more complicated to analyze. Thus 
the subject of associated production of a stop pair with a jet or Higgs requires further 
detailed studies. 

Alternative probes of light stop scenarios involve search for the heavier stop through 
its t 2 —)• tih/Z decay. Indeed the branching ratios for these channels in particular, for 
the decay t 2 —)■ tiZ, are expected to be large because of the observed large Higgs mass. 

Interestingly, for the allowed points in our light stop scenario the heavy Higgs 
phenomenology is found to be very interesting. Searches of the heavy Higgses provide 
additional channels to probe the model, either though conventional signatures in the rr 
decay channel or taking advantage of decays into electroweakinos, including invisible 
decays. Our MSSM scenarios, with no relationships between different gaugino masses, 
in fact allows for considerable values for these branching ratios and yet satisfy the LHC8 
constraint on electroweakinos. 

Finally we highlight the complementarity with DM searches, nearly all points with 
the stop-NLSP mass difference below the W mass will be tested at XenonlT provided 
the NLSP is a chargino, this means that one region that is hard to cover at the LHC 
via squark and electroweakino searches will nicely be probed by ton scale detectors. 
Unfortunately, it is much harder to cover the region with a stop NLSP with direct DM 
searches. 

In this analysis we have explicitly rejected the long-lived particles, however in our 
initial sample a large fraction of the points involved long-lived charged particles, in 
particular charginos. Existing searches for long-lived particles constrain severely the 
dominantly wino charginos in certain mass ranges [218]. Clearly improved analyses 
could provide a handle to probe this region of the parameter space left unexplored in 
this work. 

Thus, we have shown that the light stops, being actively hunted at the LHC, can 
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be probed by more than one means. While the LHC Rnn2 will bring interesting resnlts 
and hopefnlly a BSM signal, possible ways of constraining the light stop scenario via 
indirect constraints shonld also be considered. 
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